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Abstract 

We exhibit smooth initial data for the 2D water wave equation for which we prove 
that smoothness of the interface breaks down in finite time. Moreover, we show a stability 
result together with numerical evidence that there exist solutions of the 2D water wave 
equation that start from a graph, turn over and collapse in a splash singularity (self 
intersecting curve in one point) in finite time. 
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1 Introduction 

We consider the 2D water wave equation, which governs the motion of the interface between 
a 2D inviscid incompressible irrotational fluid and a vacuum, taking gravity into account 
but neglecting surface tension. We prove that an initially smooth interface may in finite 
time become singular by the mechanism illustrated in fig. [TJ We call such a singularity a 
"splash" . We also present numerical evidence for a scenario in which the interface starts out 
as a smooth graph, then "turns over" after finite time, and finally produces a splash, as in 




Figure 1: Splash singularity. The interface collapses at one point. 



The equations of motion in M 2 for the density p = p(x,t), x E M 2 , t > 0, the velocity 
v = (v , v 2 ) and the pressure p = p(x, t) are: 



p{vt + v ■ X7v) 
Pt + v • Vp 
V • v 



-V P -(0,p), 

o, 

0. 



(1) 
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Above, the acceleration due to gravity is taken equal to one for the sake of simplicity. 
The free boundary is parameterized by 



dn j (t) = {z(a,t) = (z 1 (a,t),z 2 (a,t)) 
where the regions W(t) are defined by 

r o, xGrfft) 

p{x 1 ,x 2 ,t) - | ^ x £ n 2 (t) = R 2 - Q\t). 

We assume that the fluid is irrotational, i.e. the vorticity X7 1 - ■ v = 0, in the interior of 
each domain £l J (j = 1, 2). The vorticity will be supported on the free boundary curve z(a, t) 
and it has the form 

V" 1 • v(x, t) = oj(a, t)S(x — z(a, t)), 
i.e. the vorticity is a Dirac measure defined by 



(V ± -v,r])= / Lu(a,t)r](z(a,t))da, 



with r/(x) a test function. We present results for the following geometries: 
1. open curves asymptotic to the horizontal at infinity 

lim (z(a,t) - (a,0)) = 0. 

a— >oo 



2. periodic curves in the space variable 

z(a + 2kir, t) = z{a, t) + 2/cvr(l, 0). 



3. closed contours 

z(a + 2kir, t) = z(a, t). 

However, in this paper we will only deal with the first the others are similar. 

One method to derive the equations for the evolution of v, p is to write the velocity as the 
orthogonal gradient of the stream function, take the curl and recover the velocity by inverting 
the Laplacian, i.e., we apply the Biot-Savart law. Here we use the fact that the vorticity is 
concentrated on the interface; 

v(x,t) = V ± A- 1 (V ± -v)(x,t) = ±- [ { *~ z{ ^\ {a,t)da, 

2?r J r \x - z(a,t)\ 2 

with x / z(a, t). 

Taking limits of the above equation, by approaching the boundary in the normal direction, 
we obtain the velocity of the interface, to which we can add any term c in the tangential 
direction z a without modifying the geometry of the interface. Thus the interface satisfies 

z t (a, t) = BR(z, uj)(a, t) + c(a, t)z a (a, t), (2) 
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where the Birkhoff-Rott integral is defined by 

2vr y R |z(a,t) - z(A*)l 
The system is closed by using Euler equations (for details see for example [S]); 



M 2 



uj t {a,t) = -2d t BR(z,u))(a,t) ■ d a z(a,t) - d a (— — ^)(a,t) + d a (cui)(a,t) 

+ 2c(a,t)d a BR(z,oj)(a,t) ■ d a z(a,t) + 2d a Z2{a,t). 

Then the dynamic equations, for the interface z = z(a,t) and the vorticity are the system 
given by @ and ^ and are known as the water-wave equations. 

Taking the divergence of the Euler equation and recalling that the flow is irrotational 
in the interior of the regions we find that 

-Ap = \Vv\ 2 > 

which, together with the fact that the pressure is zero on the interface implies by Hopf's 
lemma in Q 2 (t) that 

a(a,t) = -\z^(a,t)\d n p(z(a,t),t) > 0, 

where d n denotes the normal derivative. This is known as the Rayleigh- Taylor condition (see 
[20] and [22]) which was first proved by Wu in [23] and (23]. 

The first results concerning the Cauchy problem for the linearized version of water waves 
and small data in Sobolev spaces are due to Craig [11], Nalimov [IT], Beale et al. [5] and Yoshi- 
hara |27| . The well-posedness in Sobolev spaces for the water-wave problem was proven by Wu 
in [23] with the assumption that the initial free interface is non-self-intersecting (satisfies the 
arc-chord condition) . For recent work on local existence see Wu [21] , Christodoulou-Lindblad 
[S], Lindblad [TO], Coutand-Shkoller [TO], Shatah-Zeng [21j, Zhang-Zhang [2E], Cordoba et al. 
[9], Lannes [13], [TO], Alazard-Metivier [1] and Ambrose- Masmoudi [3j. The issue of long time 
existence has been treated in Alvarez-Lannes [2] where well-posedness over large time scales 
is shown and different asymptotic regimes are justified. Wu proved in |25J exponential time 
of existence in two dimensions for small initial data, and Germain et al. in |12j and Wu in 
[26j global existence for small data in the three dimensional case (two dimensional interface). 
In [6] and [7J, Castro et al. showed that there exist large initial data parameterized as a 
graph for which in finite time the interface reaches a regime in which it is no longer a graph. 
For previous numerical simulations showing this phenomenon see Beale et al. [1]. 

The outline of the paper is the following: in section 2 we will describe the equations 
in a transformed domain which will circumvent the problem of having a singularity where 
the arc-chord condition fails as the curve self-intersects, i.e. a splash singularity forms. In 
section 3 we will outline the proof of a local existence theorem for the equations in the new 
domain, both for the analytic case and for Sobolev spaces. Section 4 will be devoted to a 
stability theorem, whereas section 5 will comment on the numerical results obtained towards 
the splash singularity starting from a graph. Finally in section 6 we describe the ideas that 
we hope will lead to a computer assisted proof of the existence of a solution that starts as a 
graph and ends in a splash. 
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2 The equations in the tilde domain 



In this section we will rewrite the equations by applying a transformation from the original 
coordinates to new ones which we will denote with a tilde. The purpose of this transformation 
is to be able to deal with the failure of the arc-chord condition. We start by reformulating 
the set of equations, in the non tilde domain, for the case of a periodic contour in terms of 
the velocity potential. From (fij) and since v is irrotational in 2 (t) we have that: 



V x v = in n 2 (t) 
V-v = mQ 2 (t) 

<f>t(x,y,t) + ^\v(x,y,t)\ 2 = -(p(x,y,t) -p*(t)) -y in Sl 2 (i) 

zt(a, t) = u(a, t) + c(q, t)z a (a, t) 
$ a (a,t) = u(a,t) ■ z a (a,t) 
p = constant = in 
z(a,0) = z°(a) 

d> Q (a,0) = $°(a) (4) 

where <f> is the velocity potential, $(a,£) is its limit at the interface coming from the fluid 
region, V0 = v, p* is a function of t alone, c is a free quantity which represents the reparam- 
eterization freedom, u{a, t) is the limit of the velocity at the interface coming from the fluid 
region and &a( a ) nas zero mean - We also want the velocity to be in 1? {Vl?{f) n ([— vr, it] x 
and periodic in the x— coordinate i.e. 

v(x + 2-7T, y) = v(x, y) in Q 2 (t). 

Note that $ is periodic in the horizontal variable, because v is periodic and v(x, y) tends 
to zero as y tends to — oo. 

In order to simplify this system we use the stream function tjj and consider the equations 
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tp(x + 2n,y,t 

z t (a,t 
4>(x,y,t 

z(a, 
<S> a {a,0 



= in n 2 (t) 

\z a (a,t)\ 
= ip(x, y, t) in tt 2 (t) 
is 0(1) as y — > — oo 
= in tt 2 (t) 

= u(a, t) + c(a, t)z a (a, t) 
is the harmonic conjugate of ip(x,y,t) in Q 2 



1 



u(a, t)\ 2 + c(ct, t)u(a, t) ■ z a (a, t) — Z2(a, t) + p*(t) 



2 

z°{a) 



(5) 



Although we may take as an initial condition the tangential component of the velocity 
multiplied by the modulus of the tangent vector, i.e. $ a , we can also solve the system §5§ by 
taking as an initial condition the normal component of the initial velocity multiplied by the 
modulus of the normal vector, i.e. (^(a, t) = ip(z{a, t),t)), as we can transform one into 
the other. 

It can be checked that solutions of the system ^ are also solutions of the system Q. 
Let us consider z(a,t) = (zx(a, t), Z2(a, t)) = P(z(a,t)) where P is a conformal map defined 
in the water region that will be given as: 



P(w) 



tan 



1/2 



w G C, 



for a branch of the square root that separates the self-intersecting points of the interface. 
Here P(w) will refer to a 2 dimensional vector whose components are the real and imaginary 
parts of P(wi + W2). In this setting, P^ 1 (z) will be well defined modulo multiples of 2n. 

The water wave equations are invariant under time reversal. To obtain a solution that 
ends in a splash, we can therefore take our initial condition to be a splash, and show that there 
is a smooth solution for small times t > 0. As initial data we are interested in considering 
a curve that intersects itself at one point, as in fig. |4j More precisely, we will use as initial 
data splash curves which are defined as follows: 

Definition 2.1 We say that z{a) is a splash curve if 

1. Z\{ot) —a,Z2(a) are smooth functions and 2ir -periodic. 

2. z[a) satisfies the arc-chord condition at every point except at ct\ and ot2, with a\ < «2 
where z(ot\) = z{pi2) and |z a (ai)|, |^ a (a2)| > 0. This means z{ct\) = z{a.2), but if we 
remove either a neighborhood of ot\ or a neighborhood of a.2 in parameter space, then 
the arc-chord condition holds. 
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3. The curve z(a) separates the complex plane into two region; a connected water re- 
gion and a vacuum region. The water region contains each point x + iy for which 
y is large negative. We choose the parametrization such that the normal vector n = 

(—dc,z 2 (a),d a z 1 (a)) ■ , , ,i 

i Id — ^—^ points to the vacuum region. 

4- We can choose a branch of the function P on the water region such that the curve 
z(a) = P(z(a)) satisfies: 

(a) Z\{ot) and Ziipt) are smooth and 2ir -periodic. 

(b) z is a closed contour. 

(c) z satisfies the arc-chord condition. 

We will choose the branch of the root that produces that 

lim P(x + iy) = -e _i7r/4 

independently of x. 

5. P(w) is analytic in w and j£(w) / i/ro belongs to the water region. 

6. z{a) / q l for I = 0, ...,4, where 



,o = ( „,0), ""={t2^)- q ' = (T2'T2 

(6) 

From now on, we will always work with splash curves as initial data. Condition 6 will be 
used in the local existence theorems and can be proved to hold for short enough time as long 
as the initial condition satisfies it. We will also need that the interface passes below the points 
(±7r, 0) (or, equivalently, that those points belong to the vacuum region) in order for the tilde 
region to be a closed curve and the vacuum region to lie on the outer part. For a splash curve 
this is trivial from the definition. For more information about the transformation of both 
regions, check the figures [3] and [4] and notice that we rule out the scenarios in fig. [2} 




Figure 2: Two examples of non-splash curves. 
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We will now write ^ in the new tilde coordinates. We define the following quantities: 

ip(x, y, t) = il)(P~ l {x, y, t)), cf>(x, y, t) = ^>(P _1 (x, y, t)), v(x, y, t) = V(j)(x, y, t). 

Let us note that as ip and <p are 2ir periodic, the resulting ip and <p are well defined. We do 
not have problems with the harmonicity of ip or (j) at the point which is mapped from minus 
infinity (which belongs to the water region) by P as 4> and ifi are well defined at infinity. Also, 
the periodicity of 4> and ip causes 4> and ip to be continuous (and harmonic) at the interior of 
P(H 2 (t)). 

Let us assume that there exists a solution of ^ and that we take u n = j^r such that 
v-ufai), u n (a2) < for all < t < T, T small enough, thus z(ct, t) satisfies the arc-chord 
condition and does not touch the removed branch from P(w). 

Thus the system ([5]) in the new coordinates reads 

A$(x,y,t) = Q mP(n 2 (t)) 

^ , N = ~ l~ ( *\\ 

z(a,t) \z a {a,t)\ 

v = V ± iJj in P(n 2 (t)) 
zt(a, t) = Q 2 (a, t)u(a, t) + c(a, t)z a (a, t) 

$t(a,t) = -Q 2 (a,t)\u(a,t)\ 2 + c(a,t)u(a,t) ■ z a (a,t) - P 2 1 (S(a, t)) +p*{t) 
z(a,0) = z°(a) 

$ a (a,0) = $°(a) = $°(a). (7) 

where u is the limit of the velocity coming from the fluid region in the tilde domain and 
Q 2 ( a ,t) = \ d ^(z(a,t))\ 2 . 

We can therefore solve the Neumann problem for the stream function in the fluid domain 
with boundary conditions on the interface. In fact there exists a function Co satisfying 

i r 

ip(x,t) = — / log(|x — z{a,t)\)u(a, t)da, x belonging to the fluid region, 
2vr J_ n 

which implies 

v(x,t) = V^(x,t) = — P-V ~ u>(a,t)da. 

2tt J_ n \x - z{a,t)\ z 

Taking limits from the fluid region we obtain 

u(a,t) = BR(z,£j) + - 2 z a . 

Z \ Zfy 
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The evolution of ui is calculated in the following way. First, let us recall the equations 



z t (a 


t) 


= Q 2 (a, t)u(a, t) + c(a, t)z a (a, t) 


®t(et 


t) 


= -Q 2 (a, t)\u(a, t)\ 2 + c(q, t)u(a, t) ■ z a (a, t) — P 2 _1 (z(a, t)) 




f) 


= u(a, t) ■ z a (a, t) 


z(a, 


0) 


= z°(a) 




0) 


= *°(a) = *«(a) 



Those equations will be the ones used in section 5 with c(a, t) = 0. Substituting the 
expression for u(a,t) and performing the change c(a,t) = c(a,t) + ^Q 2 (a, t) jjrr^rfjp we 
obtain 



zt(a, t) = Q 2 {a, t)BR(z, u))(a, t) + c(a, t)z a (a, t) (9) 
and the evolution equation for u 

£j t {a,t) = -2d t BR(z,u){a,t) ■ z a (a,t) - \BR(2,u)\ 2 d a Q 2 (a,t) 

( Q 2 (a,t) u(a,t) 2 \ _ 
- <9 Q ( ^ , L 2 ) + 2c(a,t)d a BR(z,u) ■ z a (a,t) 

+ d a {c(a,t)u(a,t))-2d a {P^ l {z{a,t))) . (10) 

Remark 2.2 Equation (10) is analogous to In fact, if we set Q = 1 m (10) we recover 
<§. 

Note that for the tilde domain, the Rayleigh- Taylor condition is the same as in the first 
domain, i.e: 

Vp(a, t) ■ Za(a,t) = Vp(a, t) ■ z a (a, t) 

where p = p o P^ 1 . 

Our strategy will be the following: we will consider the evolution of the solutions in the 
tilde domain and then see that everything works fine in the original domain. 

We will have to obtain the normal velocity once given the tangential velocity, and vice- 
versa. To do this, we just have to notice that 

<l a (a, t) = u(a, t) ■ z a (a, t) = BR{z, Co) ■ z a (a, t) + ~ ^ ' \ 

From that, we can invert the equation and get ui, plug it into the following expression for ip: 

i r 

V'w = — / logf \x — z(a>) \)cj(a)da 
2tt . 
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and restrict ourselves to the interface to get ^(a, t). Taking a derivative in a we can recover 
the normal component of the velocity. An analogous reasoning can be done to get the 
tangential velocity from the normal by solving the complementary Neumann problem for <f>. 

We now note that a solution of the system ([7]) in the tilde domain gives rise to a solution 
of the system §5§ in the non-tilde domain, by inverting the map P. In fact, this will be the 



implication used in Theorem 3.1 (finding a solution in the tilde domain, and therefore in the 
non-tilde). 

Remark 2.3 A similar argument works for the other two settings (closed contour and asymp- 
totic to horizontal) by choosing an appropriate P(w) that separates the singularity. 



3 Local existence at the splash 

The main result in this section is a local existence proof for the splash singularity. To avoid 
the arc-chord condition failure, we will prove the local existence in the tilde domain. This can 
be done in two different settings, namely in the space of analytic functions and the Sobolev 
space H s . 

Theorem 3.1 Let z°(a) be a splash curve such that z\{a) — a^z^icx) G H 4 (T). Let u°(a) ■ 
(z°)- l (q) G F 4 (T) satisfying: 

1. U ■ ^) (en) < 0, L° ■ ^] (aa) < 0, 

2. I u°- ^y^ds = [ u°{a) ■ {zD^da = 0. 

J da \ z a\ Jl 

Then there exist a finite time T > 0, a time-varying curve z(a,t) 6 C([0, T]; i? 4 (T)) satisfy- 
ing: 

1. z\{a,t) — a,Z2(a,t) are 2ir-periodic, 

2. z(a,t) satisfies the arc-chord condition for all t G (0, T], 

and u(a,t) G C([0, T]; H 3 (T)) which provides a solution of the water wave equations §5§ with 
z(a,0) = z°(a) and n(a,0) • (z a )- L (a,0) = u°{a) ■ O^-^a). 

Sketch of the proof: Using the fact that there is local existence to the initial data in the tilde 
domain and applying P^ 1 to the solution obtained there, we can get a curve z(a,t) that 
solves the water waves equation in the non tilde domain. This leads to the proof of Theorem 
13. li Details on the local existence in the tilde domain are shown below. 

3.1 Local existence for analytic initial data in the tilde domain 

In this subsection, we will work on the tilde domain, and all tildes will be dropped for the 
sake of simplicity. 
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We will work with c = 0. We have the following system: 

2 



I dP_ 
I dw 



Zt 

u 

2 



{p-\z(a,t))) 
P 2 l (z{a,t)) 



l\ d £(p-\z))\ 2 M 2 - 

BR(z,u>) + ^z a 
| + BR(z,uj) ■ z a 

J_ l+(z 1 (a,t)+iz 2 (a,t)) 4 
16 Zi(a,t)+iz2(a,t) 
i+(z 1 (a,t)+iz2(a,t)) 2 
i—(z\(a,t)+iZ2{a,t)) 2 



(11) 



In 



We demand that z°(a) / (0,0) to find the function ( ^ } {P~ 1 (z{a,t))) well defined. This 
condition is going to remain true for short time. We also consider z°(a) ^ q\ I = 1, 4 in 
Q to get P% (z(a,t)) well defined. Again this is going to remain true for short time. 

We consider the space 



H 3 (dS r ) = {/ analytic in S r = {a + irj, \rj\ < r} : 



. < CO 



(asv) + \Kf\\ 



where 



2 

L 2 (dS r ) 



/ I /(a db ir)\ 2 da, f 27r-periodic 



and (z, $) G (H 3 (dS r )) 3 = X r . 
We have the following theorem: 



,0. Zo 



a 



Theorem 3.2 Zei z°(a) 6e a splash curve and let u -t^t\ 
velocity such that 

(z?(a) -a,^,(a),$°(a)) G X ro , 

/or some ro > 0, and satisfying: 



a) 6e i/ie initial tangential 



u 



(*£)- 



(4) 



(ai) < 0, I u 



(4)" 



(«2) < 0, 



u°(a) • A* = 0- 



'dfi Pal 

Taen i/tere exisi a finite time T > 0, < r < ro, a time-varying curve z(a,t) satisfying: 

1. P~ l {z\{a,t)) — a, P~ 1 (z2(a, t)) are 2it -periodic, 

2. P~ 1 (z(a,t)) satisfies the arc-chord condition for all t G (0,T] ; 
and n(a, t) wi/i 

(zi(a,t)-a,^(a,t),$ tt (a,t)(a,t)) G C([0,T],X r ) 

which provides a solution of the water waves equations (11) with z°(a) = P(z°(a)) and 
u(a,0) ■ (z a ) ± (a,0) = u°(a) ■ (*°)£(a). 



The main tool in the proof is to use an abstract Cauchy-Kowaleswki theorem from [18 
and |19j as in for example [6]. 
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3.2 Local existence for initial data in Sobolev spaces in the tilde domain 

We will take the following c(a,t): 



We will also define an auxiliary function <p(a, t) analogous to the one introduced in [5j 
(for the linear case) and [3] (nonlinear case) which helps us to bound several of the terms 
that appear: 

Q 2 {a, t)u}(a, t) . , 
<p{a,t ) = * V ' / \ - - c(a,t)\z a (a,t)\. 12 
2\z a (a,t)\ 

Theorem 3.3 Let z°(a) be a splash curve such that z\(a) — a, z^ipi) E H 4 (T). Let u°(a) ■ 
(zl)^{a) E F 4 (T) satisfying: 

1. U ■ ^) (ai) < 0, L° ■ ^) (as) < 0, 

2. f u°- ^Lds = [ u°(a) ■ {zD^da = 0. 

JdCl \ z a\ Jt 

Then there exist a finite time T > 0, a time-varying curve z(a,t) E C([0,T]; H 4 ) satisfying: 

1. P^ 1 (zi(a, t)) — a, P~ 1 (z2(a, t)) are 2ir -periodic, 

2. P~ 1 (z(a,t)) satisfies the arc-chord condition for all t E (0,T], 

and u(a,t) E C([0, T]; H 3 (T)) which provides a solution of the water waves equations ^ 
with z°(a) = P(z°(a)) and u(a, 0) • (5 a ) ± (a, 0) = u°(a) • (5°)£(a). 

Sketch of the Proof: 

In the proof, for the sake of simplicity, we will drop the tildes from the notation. 

The proof will use the properties of c(a,t) and (p(a,t) to get an extra cancellation to 
help us derive energy estimates. Moreover, this choice of c will ensure that the length of the 
tangent vector of z(a,t) depends only on time. 

Here we define the energy E(t) by 

E(t) =\\4h*V) + I ^Kz\ 2 da +\\F(z)f LOB (t) 

Jt w\ 

| 2 12 4 1 



n 



where the L°° norm of the function 

F(z) = 



z(a, t) — z(a — /3,t)\ 



measures the arc-chord condition, 



a z = [ BR t {z,u) + -^-BR a (z,ojU . ^ + 

I ~-n I / ^ -^a | 



Zat + i — r£ 



+ g 



BR(z,uj) + 



a; 



2 1 £rv I 2 



(13) 



(VQ)(z)-4 + (VP 2 - 1 )(z)-^ 



is the Rayleigh- Taylor function, 



and finally 



m (Q (T z )(t) = mmQ (a,t)a z (a,t), 



m(q l )(t) = min |z(a, i) — q l \ 



for Z = 0, ...,4. We proceed as in [9]: The bound for the operator ( I + J) , where Jo; = 
2BR(z,uj) ■ z a , and some rather routine estimates allow us to find 

+ \\F(z)\\Ut) + MUt) + h\\Ut) + wtS!^ + E ^7^) * 



for C and A; universal constants. Above we use that 

l2„ U |2 U |2 



Further we obtain 
where 



*(J z \Z a \ I (-,4 1 2 7 ^ 

\otz\ da. < 



Q 2 cr z \z a \ 2 ' a 1 - m(QV. 
d_ / /" QV 



I ^\dU 2 da<E 2 (t). 

JT \ z a\ 



(J T j^\dU 2 da)<CE k (t) + S(t) 
S(t) = 2 / ^a^z • iff(^)da. 



We use (12), (10) and (13) to get 



Q 2 <Jz 



where "control" is given by lower order terms and unbounded terms (such us d^ip) that can 
be estimated with energy methods in terms of E(t). Therefore, it allows us to get 

d 

~dV 

which together with above inequalities yields 

d 



\^ 1/2 d 3 M\h(t) < CE k (t) - S(t) 



dt 



E(t) 



< CE k (t). 



Local existence follows using standard arguments with the apriori energy estimate. 



□ 
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4 Structural Stability 



Again, in this section, we will omit the tildes from the notation. This section is devoted to 
establish a stability result. It will allow us to conclude the following: if (x, 7) approximately 
satisfies equation (14) , then near to (x, 7) there exists an exact solution (z,uj). Below is the 
theorem. 

Theorem 4.1 Let 

D(a, t) = z(a, t) — x(a, t), d(a, t) = w(a, t) — 7(0;, t), TJ(a, t) = ip(a, t) — ((a, t) 
where (x, 7, £) are the solutions of 



b 



Q 2 {x)BR(x,-f)+bx a + f 
a + 7r 



2vr 



" (Q 2 BR(x, j)) a -^-da - f (Q 2 BR(x, 7 )) /3 ^L_d/3 

TV \^a\ J— 7T \^ot\ 



+ 



a + it 
2tt 



n JtJ \ x p\ 



it 



+2BR t (x, 7) • x a = -{Q 2 {x)) a \BR(x, 7 )| 2 + 2bBR a (x, 7) • x a + (6 7 ) Q 
-2{P£ x {x)) a + g 
C(«.*) = ^M a -M°,*)Ma,t)|, 



(14) 

(2, w, (/?) are the solutions of ( 14 ) with f = g = and £ i/ie following norm for the difference 



Then we have that 



<C(t)(£(t) + S(t)) 



where 



and 



Cit) = C(E(t), \\x\\ H&+h (t), hH^x (t), HCll^i (t), \\F(x)\\ L ~ (t)) 



m = (\\f\\„ 5+ l(t) + \\9\\„ 3+ l (*))* + 



i(*) + ll5ll^3 + i(i)) 2 , 



6z<7 enough. Here E(t) is defined in the proof of Theorem 3.3. 

Sketch of the Proof: The equation for (x, 7, £) is the same as the one for (z,ui,(p) but for / 
and <?. The function 6 is chosen in such a way that |ac a | only depends on time. Then it allows 
us to get the following estimates: 



d_ 
Jt 



\DVn* + 



+ ll^l 



L 2 , 



<C(t)(£(t) + S(t)). 
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Further we obtain 



with 



For T> one finds that 



( jf ^|^£>| W) < C(t)(£(t) + 6(t)) + S(t) 
S(t) = 2 f ^d^D ■ ^H{d A a V)da. 

JT \ z a\ \Z a \ 



dlV t = -?°*rdiD ■ i + "control" 



i 1 2 or | 



where "control" denotes terms which can be estimated by C(t)(£ (t) + S(t)), which yields 

|||AV^3p||2 2(t) < C (t)(£(t) + 5(t)) - S(t). 
Then the desired estimate follows. □ 



5 Numerical results 

In order to illustrate the splash singularity, several numerical simulations were performed. 
The simulations were done following the scheme proposed by Beale, Hou and Lowengrub 
[I] adapted to the equations on the tilde domain (i.e. taking into account the impact of Q 
on the equation). Instead of having an evolution equation for u, they introduce a velocity 
potential <ft and study its evolution through time subject to the constraint imposed by being 
a potential. This is the set of equations Q. The initial data on the non-tilde domain was 
given by: 

*?(a) = a + ^ (- y - 1.9^ sin(a) + ^ sin(2a) + i (| - 1.9) sin(3a) 

13 1 

zl{a) = — cos(a) - — cos(2a) + — cos(3a) 

Note that z (|) = z (— f ) (splash). Instead of prescribing an initial condition for uj, we 
prescribed the normal component of the velocity to ensure a more controlled direction of 
the fluid. From that we got the initial w(a,0) using the following relations. Let ip be such 
that V^tp = v and ^(q) its restriction to the interface. Recall that we can transform the 
initial condition on the normal component of the velocity into an initial condition on the 
tangential component by applying the transformations described in section 2. The initial 
normal velocity is then prescribed by setting 

u^(a)\z a (a)\ = ^a{a) = 3 • cos(a) — 3.4 • cos(2a) + cos(3a) + 0.2cos(4a). 

The simulations were done using a spatial mesh of N = 2048 nodes and a time step 
At = 10 -7 . The time direction was set to run backwards (from the splash to the graph) and 
the graph was obtained at approximately T g = 6.5 • 1CP 3 . Note that the normal component 
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Figure 3: Tilde domain at times t = (Red - splash. The splash points are located on 
the open segment joining the points (v2/2, v2/2) and (— v2/2, — v2/2), approximately at 
(±0.27284156, ±0.27284156)), t = 4 ■ 10" 3 (Blue - turning) and t = 7 ■ 10~ 3 (Black - graph) 




Figure 4: Zoom of the splash singularity at times t = (Red - splash), t = 4 • 10 3 (Blue - 
turning) and f = 7 • 10~ 3 (Black - graph) 



3.1 



as 



of the velocity (±f ) > at the splash, which satisfies the hypotheses of Theorem 
we are running time backwards. Getting the potential of the initial condition from loq and 
zq and the transformation of all the initial data to the tilde domain is a trivial computation 
taking care to choose the appropriate branch of the square root. See figures [3] and [4} 



6 Further research 

We would like to exhibit a water-wave solution whose interface starts as an i? 4 -smooth graph 
at time zero, and ends in a splash at time T. We sketch a few ideas that may lead to a 
rigorous computer-assisted proof of the existence of such a solution. We will work in the tilde 
domain. 

A simulation as in Section 5 leads to an approximate solution x(a,t), j(a,t) with the 
desired properties. Thus x(-,t) describes a graph when t = and a splash when t = T. 



Moreover, we believe that equations similar to ( 14 ) hold, with very small / and g 



We may suppose that x and 7 are known piecewise-polynomial functions on [0, 2ir] x 
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[0, T]. Using interval arithmetic [16], one can compute rigorous upper bounds for appropriate 
Sobolev norms of / and g. We hope that these upper bounds will be very small. 

Next, we solve the water-wave equations for z,u, starting at time T, and proceed- 
ing backwards in time. We take our initial (z,uj) at time T to be a splash, very close to 
T), 7(-, T)) in a high Sobolev norm, [j] 

We want to compare the exact solution (z,uj) with the approximate solution (x,7), using 



the quantity £(t) as in Theorem 4.1 Since (z,ui) and (a;, 7) are very close at time T, we will 
be able to show easily that 

£(T) < £1, for a very small, computable constant E\. (15) 

Moreover, the functions x and 7 are known; and we also know upper bounds for Sobolev 
norms of / and g. Therefore, the ideas in the proof of Theorem |4.1[ together with interval 
arithmetic, should lead to a rigorous proof of the differential inequality 



<C,£(i)+£ 2 , (16) 



where C\ and £2 are computable constants. 

We hope that £2 is very small, since / and g have small Sobolev norms; and we hope that 
C\ won't be too big. 



Once we establish (15) and (16), we will then know that our water wave solution (z,uj) 
exists for all time t S [0, T], and that £(0) < £3 for a (hopefully small) computable constant 
£3- 

From the definition of £ (t), we will then easily deduce that z(-,0) — x(-,0) has norm at 
most £4 in H 4 (M./2it'Ij), for a computable constant £4. 

If £4 is small enough, this in turn implies that the interface z(-,0) is an ff 4 -smooth graph. 
Thus {z,uj) is an exact solution of the water-wave equation, whose interface is an if 4 -smooth 
graph at time 0, and a splash at time T. 

We hope that a proof along these lines can be made to work. 
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